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We  are  studying  a  subwavelength  slit  in  a  lossy  metal  film  illuminated  by  electromagnetic  waves.  Periodic 
grooves  are  placed  on  the  output  side  surrounding  the  slit.  Surface  plasmons  along  the  metallic  surface  on  the 
output  side  are  excited  by  the  electromagnetic  wave  propagating  through  the  slit.  The  grooves  modulate  the 
surface  plasmons  and  couple  them  into  propagation  modes,  which  facilitates  the  beaming  of  the  output  field. 

We  developed  an  efficient  method  to  determine  the  geometric  parameters  of  the  grooves  that  are  necessary  to 
achieve  oblique  beaming  at  any  angle  between  0  and  70  degrees.  We  also  designed  a  frequency  splitter  by 
setting  appropriate  parameters  of  the  grooves. 
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Normally  the  light  coming  out  a  subwavelength  wave¬ 
guide  will  be  diffracted  to  all  angles.  With  the  help  of  surface 
waves  along  the  output  surface,  the  emitted  light  can  be  con¬ 
fined  in  a  small  angle  range.  This  phenomenon  is  called 
beaming  and  was  discovered  in  2002. 1  Since  then,  it  has 
been  studied  extensively  experimentally  and  theoretically  us¬ 
ing  metallic  structures 1-6  or  photonic  crystals. 7-11 

Beside  the  research  on  forward  beaming  mentioned 
above,  oblique  beaming  was  also  demonstrated  using  either 
metallic  structures 12-14  or  photonic  crystals.15  Obviously  it  is 
important  to  control  the  beaming  angle  of  the  oblique  beam¬ 
ing  freely.  In  this  paper,  we  present  an  efficient  method  to 
design  a  metallic  structure  which  will  steer  the  emitted  field 
to  any  given  angle.  We  also  find  the  best  beaming  happens  at 
an  oblique  direction. 

Figure  1  shows  the  structure  studied  in  the  paper:  a  me¬ 
tallic  subwavelength  slit  surrounded  by  grooves.  It  is  an  im¬ 
portant  plasmonic  structure. 16  As  a  two-dimensional  struc¬ 
ture,  it  can  be  solved  analytically  and  serves  as  a  theoretical 
model  to  understand  the  physics  of  beaming.2  It  is  a  multi¬ 
functional  beaming  structure.  It  can  achieve  intermediate- 
field  beaming,3  far-held  beaming,2  and  specially  far-held  ob¬ 
lique  beaming  which  will  be  studied  here.  The  near-held 
focusing  was  also  implemented  by  similar  structures.20 
Grooves  at  the  input  surface  can  enhance  the  transmission 
dramatically,17,18  so  this  structure  can  improve  transmission 
and  directionality  of  the  output  held  simultaneously.12  Re¬ 
cently  this  structure  was  integrated  with  a  semiconductor  la¬ 
ser  to  converge  the  emitted  fields.19 

Since  electromagnetic  helds  decay  rapidly  inside  the 
metal,  the  metal  can  be  “replaced”  by  surface  impedance 
boundary  conditions  (SIBC),21  which  states  the  electromag¬ 
netic  helds  at  the  metal  surface  satisfy 


E\\  =  Zn  X  Hh  Z  =  \/imp.0/eme0. 

Here  n  is  the  normal  unit  vector  pointing  out  of  the  metal; 
£11  and  //||  are  components  of  electric  and  magnetic  helds 
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FIG.  1.  (Color  online)  Schematic  representation  of  the  subwave- 
length  metal  slit  surrounded  by  upper  and  lower  periodic  grooves 
embedded  in  free  space. 
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FIG.  2.  (Color  online)  (a)  and 
(b)  Amplitudes  and  phases  of 
fields  at  symmetric  grooves  with 
a=b  =  500  nm;  (c)  and  (d)  the  lin¬ 
ear  fitting  results  of  phases  of  the 
fields  at  grooves  for  different  a 
when  (c)  ft  =  500  nm  and  for  dif¬ 
ferent  b  when  (d)  a  =  500  nm. 


parallel  to  the  metal  surface.  Z  is  the  impedance  of  the  metal; 

and  jjim  are  relative  permittivity  and  permeability  of  the 
metal;  e0  and  //„  are  vacuum  permittivity  and  permeability. 

The  coordinate  system  is  defined  in  Fig.  1.  x=0  is  the 
output  surface  of  the  slit.  When  the  incident  wave  is  TE 
polarized  (the  nonzero  components  of  the  electromagnetic 
fields  are  Hz,  Ex,  and  Ey),  surface  plasmons  exist  at  the  metal 
surface.  In  this  paper,  we  focus  on  the  fields  at  and  after  the 
output  surface  of  the  slit.  At  the  output  surface,  the  SIBC 
becomes 


Ey  +  ZHz  =  0. 

Define  F=Ey+ZH~.  Given  the  F  field  along  x=0,  we  can 
calculate  the  electric  and  magnetic  fields  at  any  points  in  the 
x>0  half-plane.  Define  the  far-held  angular  transmission 
P(9)  as 

P{6)  =  lim  r|S(r)|,  (1) 


where  S{r )  is  the  Poynting  vector  at  the  point  r 
=  (r  sin  0,r  cos  6).  Then  the  total  transmission  T  is 


T  = 


p{e)dd , 


and  P{0)  can  be  calculated  from  the  F  fields  as 


P{  0)  =  — 
2Z0 


cos  6 


cos  9  +  Z/Z, 


2  p® 

dy'F{0,y')e-b 

J  —00 


ik o  sin  0y' 


(2) 


Here  we  assume  the  x>0  half-plane  is  free  space;  k0  and 
Z0  are  the  wave  vector  and  the  impedance  in  free  space. 
Then  we  define  the  normalized  far-field  angular  transmission 
as  I{6)  =  P{9) / T  and  use  it  to  describe  the  far-field  distribu¬ 
tion. 


Because  of  the  surface  impedance  boundary  condition,  the 
F  field  along  x=0  is  zero  except  the  opening  area  of  the  slit 
and  grooves.  Since  the  widths  of  the  indentations  (slit  and 
grooves)  are  much  smaller  than  the  incident  wavelength,  we 
assume  the  field  is  constant  inside  an  indentation  and  Eq.  (2) 
is  approximately  the  summation  of  several  point  sources.  We 
will  define  the  F  field  at  the  ith  grooves  as  Ft\  positive  i 
means  the  groove  above  the  slit,  negative  i  means  the  groove 
below  the  slit  and  i  =  0  represents  the  slit.  In  this  paper,  Fi  are 
obtained  by  using  the  finite  element  method  software  COM- 
SOL  Multiphysics.  They  can  also  be  calculated  by  using  a  set 
of  self-consistent  linear  equations.22-23 

By  changing  the  geometric  parameters  of  the  grooves,  we 
can  change  the  F  fields  inside  the  grooves  and  control  the 
output  field  distribution  furthermore.  It  is  impossible  to  study 
all  possible  choices  of  the  grooves.  We  focus  on  the  periodic 
grooves  in  this  paper.  As  Fig.  1  shows,  all  the  grooves  have 
the  same  shape  with  width  Wg- 40  nm  and  depth  h 
=  100  nm.  The  slit  width  is  Wx = 40  nm.  All  the  grooves 
above  the  slit  consist  in  the  upper  periodic  grooves  charac¬ 
terized  by  the  initial  location  a1  and  the  period  b\.  The 
grooves  below  the  slit  are  also  periodic  and  the  correspond¬ 
ing  parameters  are  a2  and  b2.  When  ax=a2  and  bx=b2,  the 
grooves  are  symmetric  and  we  will  use  a  and  b  to  replace  ax, 
a2  and  bx,  b2.  The  incident  wavelength  is  X  =  560  nm.  In  this 
paper,  we  try  to  find  the  best  far-field  oblique  beaming  at  any 
given  angle  by  modifying  ax,  bx,  a2,  and  b2.  Under  the  re¬ 
striction  of  periodic  grooves,  we  find  an  efficient  algorithm 
to  decide  the  parameters  with  an  explicit  physical  explana¬ 
tion. 

We  begin  our  research  by  studying  the  F  fields  when  the 
grooves  are  symmetric.  Figures  2(a)  and  2(b)  show  the  am¬ 
plitudes  and  phases  of  F,  when  a  =  b  =  500  nm.  All  the  F 
fields  are  normalized  to  the  F  field  at  the  slit  exit,  so  /-),=  1 . 
Figure  2(a)  shows  the  field  at  the  slit  is  much  stronger  than 
the  fields  at  the  grooves.  The  fields  at  the  grooves  decay 
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FIG.  3.  (Color  online)  The  fitting  error  of  the  phases  of  F,  when 
a =500  nm. 

smoothly  with  the  distance  to  the  slit  and  become  very  weak 
in  the  end.  There  are  56  grooves  in  each  side  of  the  slit  in 
this  simulation.  Further  simulations  show  adding  more 
grooves  will  not  influence  the  far-field  distribution  much. 

The  most  important  property  about  the  F  fields  is  the  per¬ 
fect  straight  lines  formed  by  the  phases  shown  in  Fig.  2(b). 
So  we  can  use  linear  functions  to  fit  them.  For  a  positive 
integer  i,  the  location  of  the  zth  groove  is  y ,  =  a  +  {i-  1  )b  and 
the  phase  is  <f>,  =  Arg(F,),  then  the  linear  function  is 

(fi,  =  kyt  +  <p0.  (3) 

This  linear  relation  can  be  explained  by  the  “single  scat¬ 
tering  approximation.”  The  surface  plasmons  are  excited  at 
the  exit  of  the  slit  and  propagate  along  the  output  surface. 
Then  they  will  interact  with  the  grooves  they  meet.  The  weak 
F  fields  at  the  grooves  suggest  the  interaction  is  weak.  Most 
of  the  surface  plasmons  transmits  through  the  grooves,  so  the 
F  fields  at  grooves  decay  smoothly.  Some  is  scattered  into 
the  free  space  and  the  rest  is  reflected.  The  single  scattering 


approximation  neglects  the  reflected  surface  plasmons.  The 
surface  plasmons  repeat  their  behavior  when  they  meet  the 
next  groove,  so  the  phase  change  between  two  neighboring 
grooves  is  constant. 

The  single  scattering  approximation  requires  that  the  fit¬ 
ting  parameters  in  Eq.  (3)  will  not  change  with  a.  Our  simu¬ 
lations  confirm  it.  The  fitting  results  are  shown  in  Fig.  2(c) 
when  b  =  500  nm  and  a  changes  from  200  to  1500  nm.  We 
can  see  clearly  that  k  is  almost  constant  and  <p0  oscillates 
weakly. 

Next  we  fix  a  =  500  nm  and  change  b  from  100  nm  to 
1500  nm.  The  slope  k  in  Eq.  (3)  is  shown  in  Fig.  2(d).  The 
wave  vector  of  the  surface  plasmons  is  ksp=l.01k0.  (The 
metal  is  gold  and  its  permittivity  is  taken  from  Ref.  24.)  We 
can  find  k  is  close  to  but  not  equals  ksp.  It  is  because  of  the 
small  phase  change  introduced  by  the  grooves.  When  b  in¬ 
creases,  which  means  less  grooves  at  the  surface,  k  con¬ 
verges  to  ksp. 

The  slope  curve  in  Fig.  2(d)  is  bumpy  when  kb~  mi  (, n  is 
an  integer),  which  is  because  of  “collective  surface  modes.”2 
When  kb  =Tm,  the  reflected  surface  plasmons  from  different 
grooves  are  in  phase,  so  the  single  scattering  approximation 
becomes  weak.  In  Fig.  3,  we  plot  the  fitting  error  of  the 
phases  when  b= 530  nm  and  400  nm.  It  is  clear  that  the 
collective  surface  mode  (£  =  530  nm)  gives  a  bigger  error. 

Now  we  have  understood  the  phase  difference  between 
two  neighboring  grooves  is  kb.  If 

kb  -  k0b  sin  0=2mi,  (4) 

where  n  is  an  integer,  all  the  grooves  above  the  slit  have 
constructive  far-field  interference  at  the  angle  9.  Figure  4(b) 
shows  the  constructive  interference  angles  of  the  upper  peri¬ 
odic  grooves  at  different  periods  calculated  from  Eq.  (4) 
based  on  the  k  in  Fig.  2(d). 


FIG.  4.  (Color  online)  (a)  Nor¬ 
malized  far-field  angular  transmis¬ 
sion  when  a=b= 500  nm;  (b) 
Constructive  interference  angles 
of  the  upper  periodic  grooves  at 
different  groove  periods  calcu¬ 
lated  from  Eq.  (4)  when  n 
=  1 , 2 ,  •  •  • ,  5  ;  (c)  and  (d)  Normal¬ 
ized  far-field  angular  transmission 
when  £  =  446  nm  and  £  =  710  nm. 


045427-3 


W.  DAI  AND  C.  M.  SOUKOULIS 


PHYSICAL  REVIEW  B  82,  045427  (2010) 


FIG.  5.  (Color  online)  (a)  and 
(b)  Amplitudes  and  phases  of 
fields  at  grooves  above  the  slit 
when  the  grooves  are  symmetric 
with  o1  =  a2  =  920  nm,  bl=b2 
=446  nm  and  when  there  are  only 
upper  periodic  grooves  with  ci\ 
=  920  nm  and  7>1  =  446  nm;  (c) 
and  (d)  Normalized  far-held  angu¬ 
lar  transmission  generated  by  the 
slit  with  upper  periodic  grooves 
alone  when  b  j=446  nm  and  b\ 
=  710  nm. 


We  also  calculate  1(0)  from  Eq.  (2)  when  a=b= 500  nm 
and  the  results  are  shown  in  Fig.  4(a).  The  normalized  far- 
held  angular  transmission  distribution  curve  has  two  peaks  at 
±6°,  which  equal  the  constructive  interference  angles  of  the 
upper  and  lower  periodic  grooves.  When  b  =  530  nm,  the 
fitting  error  is  bigger  (shown  in  Fig.  3).  But  Eq.  (4)  still  gives 
the  right  beaming  angle,  which  is  around  0°. 

Besides  the  two  peaks,  the  curve  in  Fig.  4(a)  also  has  a 
smooth  base,  which  is  the  contribution  of  the  slit.  The  exis¬ 
tence  of  the  base  means  the  energy  is  diffracted  to  all  angles. 
So  the  strong  hied  in  the  slit  deteriorates  the  directional 
emission.  However,  it  also  bring  us  the  opportunity  to  control 
the  peak  heights.  Considering  the  far-held  interference  be¬ 
tween  the  slit  and  the  upper  periodic  grooves  at  angle  0,  the 
phase  difference  is  <fi0+ka-k0a  sin  0.  So  the  peak  heights 
change  with  a ,  as  shown  in  Figs.  4(c)  and  4(d). 

Figure  4(c)  also  demonstrates  the  beaming  angles  are  in¬ 
dependent  from  a.  A  beaming  angle  is  always  a  constructive 
interference  angle  of  the  upper  or  lower  periodic  grooves, 
which  is  decided  by  the  period  only.  In  Fig.  4(d),  b 
=  710  nm  and  the  upper  or  lower  periodic  grooves  have  two 
constructive  interference  angles  and  the  value  of  a  decides 
which  one  becomes  the  beaming  angle. 

Until  now  we  always  use  symmetric  grooves,  so  the  far- 
held  distributions  are  also  symmetric  with  two  peaks.  The 
ideal  oblique  beaming  has  only  one  peak.  A  nature  idea  is  to 
remove  all  the  grooves  on  the  one  side  of  the  slit.12 

Figures  5(a)  and  5(b)  show  the  amplitudes  and  phases 
of  the  F  helds  at  the  grooves  above  the  slit  when  the  slit  is 
accompanied  by  the  upper  grooves  only  or  by  symmetric 
upper  and  lower  grooves.  The  F  helds  under  the  two 
groove  settings  are  quite  close  to  each  other,  specially  the 
phases.  The  similarity  shows  the  independence  between  the 
upper  and  lower  periodic  grooves.  The  independence  is  also 
a  nature  result  of  the  single  scattering  approximation.  Since 


we  neglect  the  reflected  surface  plasmons,  the  upward  plas- 
mons  have  no  chance  to  reach  the  grooves  below  the  slit,  and 
vice  versa.  Ref.  12  also  verihed  the  independence  between 
the  upper  and  lower  grooves.  Figures  5(c)  and  5(d)  show  1(0) 
of  a  slit  with  upper  periodic  grooves  alone  when  bx 
=446  nm  and  b2  =  l 10  nm.  Comparing  with  Figs.  4(c)  and 
4(d),  we  can  hnd  that  the  absence  of  the  lower  grooves  will 
not  influence  the  peaks  generated  by  the  upper  periodic 
grooves. 

As  Fig.  4(b)  shows,  when  b j  changes  form  0.5A.  to  X,  the 
constructive  interference  angle  changes  from  -90°  to  0°.  So 
we  can  always  design  an  oblique  beaming  structure  with  the 
beaming  angle  at  any  angle.  The  structure  is  a  subwavelength 
metallic  slit  with  only  the  upper  or  lower  periodic  grooves. 
First  we  check  Fig.  4(b)  to  get  the  period  of  grooves  /?, 
smaller  than  540  nm.  Then  we  fix  the  period  and  scan  the 
initial  location  of  the  grooves  a ,  to  maximize  the  transmis¬ 
sion  peak  height. 

Figure  4(b)  also  shows  when  h\  >540  nm,  the  upper  pe¬ 
riodic  grooves  have  several  constructive  interference  angles 
and  some  of  them  are  positive.  Then  the  lower  periodic 
grooves  with  period  b2>  540  nm  can  have  negative  con¬ 
structive  interference  angles.  So  it  is  possible  that  the  upper 
and  lower  periodic  grooves  have  the  same  beaming  angles. 
Then  we  can  design  an  oblique  beaming  structure  using  both 
upper  and  lower  periodic  grooves.  Given  an  angle  0,  we 
check  Fig.  4(b)  to  find  bl  with  the  constructive  interference 
angle  0  and  b2  with  the  constructive  interference  angle  -0. 
Then  we  study  the  angular  transmission  curves  of  a  slit  with 
upper  periodic  grooves  alone.  Fix  the  period  as  bt  and  scan 
a  |  to  maximize  the  peak  at  0.  Repeat  the  process  for  b2  to 
find  the  optimal  a2.  Then  all  the  four  parameters  ax,  a2,  bl, 
and  b2  are  obtained. 

Figure  6  is  an  example  of  the  algorithm.  The  beaming 
angle  is  -14°.  From  Fig.  4(b),  we  find  the  upper  periodic 
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FIG.  6.  (Color  online)  (a)  Normalized  far- 
field  angular  transmission  of  a  slit  with  upper 
periodic  grooves  alone  when  <^  =  980  nm, 
b  i=446  nm  and  a  slit  with  lower  periodic 
grooves  alone  with  a2  =  860  nm,  b2  =  l  10  nm; 
(b)  Normalized  far-field  angular  transmission  of 
a  slit  with  asymmetric  upper  and  lower  periodic 
grooves  when  a [=980  nm,  b\  =  446  nm, 
a2  =  860  nm,  and  £>2  =  710  nm. 


grooves  with  £>[  =  446  nm  and  the  lower  periodic  grooves 
with  b2=110  nm  both  have  the  constructive  interference 
angle  at  -14°.  Then  we  fix  b ,  =  446  nm  and  calculate  the 
normalized  far- field  angular  transmission  for  different  ax. 
The  highest  peak  is  achieved  when  a,  =980  nm  and  the  peak 
is  shown  Fig.  6(a).  Repeat  the  process,  when  a2=860  nm, 
the  lower  grooves  with  £>2=7 10  nm  generate  the  highest 
peak  at  -14°.  Then  we  put  the  upper  and  lower  periodic 
grooves  together.  The  far-held  distribution  of  the  slit  sur¬ 
rounded  by  the  asymmetric  upper  and  lower  grooves  is 
shown  in  Fig.  6(b).  Now  all  the  three  components,  the  upper 
periodic  grooves,  the  lower  periodic  grooves  and  the  slit, 
interfere  constructively.  So  we  get  a  higher  far-held  transmis¬ 
sion  peak  at  0=  — 14°  than  the  peaks  in  Fig.  6(a). 

Figure  7  shows  the  properties  of  the  highest  transmission 
peaks  we  can  achieve  at  different  beaming  angles  using  up¬ 
per  periodic  grooves  only  or  using  asymmetric  upper  and 
lower  periodic  grooves  based  on  the  methods  mentioned 
above.  Adding  grooves  below  the  slit  improves  the  transmis¬ 
sion  peaks.  The  highest  peaks  are  achieved  round  20°.  When 
the  beaming  angle  is  big,  cos  0  in  the  numerator  of  Eq.  (2) 
suppresses  the  peaks.  When  the  beaming  angle  is  small,  the 
system  is  close  to  the  collective  surface  modes.  The  helds  in 
the  grooves  decay  quickly2  and  the  grooves  tend  to  cancel 
each  other  because  of  the  phases  deviation  from  the  linear 
relation.  So  the  peaks  of  the  normalized  far-held  angular 
transmission  are  low. 

The  shape  of  the  grooves  will  also  influence  the  quality  of 
the  beaming.2  The  helds  at  the  grooves  can  be  very  weak 
given  some  special  groove  shapes.  So  the  output  held  is  de¬ 
cided  by  F0  mainly  and  the  beaming  is  suppressed.  But  the 
single  scattering  approximation  tell  us  the  shape  effect  is 
independent  from  a  and  b.  So  we  choose  h  and  Wg  which 


give  fairly  strong  helds  at  the  grooves  when  a  =  b  =  500  nm. 
This  groove  shape  also  gives  us  good  beaming  for  other  a 
and  b  values  as  predicted. 

We  can  also  design  a  frequency  splitter  using  asymmetric 
grooves:  two  close  wavelengths  enter  the  slit  and  leave  at 
very  different  angles.  We  have  understood  the  upper  periodic 
groove  has  only  one  possible  beaming  angle  when  its  period 
is  smaller  than  the  wavelength  and  the  angle  satishes 

kb  -  k0b  sin  0=  2tt.  (5) 

The  phase  difference  between  the  upper  periodic  grooves  and 
the  slit  at  the  beaming  angle  is 

ka  +  </>0  -  k0a  sin  0=  liralb  +  <p0.  (6) 

We  can  hnd  two  wavelengths  Xj  and  X2  satisfying 
fA0(X|)- (/<o( X2)=  it.  When  the  upper  grooves  and  slit  inter¬ 
feres  constructively  at  one  wavelength,  they  have  destructive 
interference  at  the  other  wavelength. 

In  Fig.  8(a),  the  two  wavelengths  are  650  and  695  nm. 
When  the  depth  of  grooves  is  /t=  100  nm,  the  difference  be¬ 
tween  4>q  at  the  two  wavelengths  is  1.7.  To  get  bigger  differ¬ 
ence,  we  have  to  increase  h  to  500  nm.  The  period  of  the 
upper  periodic  grooves  is  £>[  =  500  nm.  The  constructive  in¬ 
terference  angles  of  the  upper  grooves  are  around  -17°  for 
both  wavelengths.  But  because  of  the  interference  between 
the  upper  grooves  and  slit,  only  the  field  with  \  =  650  nm 
has  a  transmission  peak  at  its  constructive  interference  angle. 
For  the  lower  periodic  grooves,  we  choose  b2  =  bx  and  a  dif¬ 
ferent  a2  to  enhance  the  peak  \  =  695  nm.  Though  the  two 
incident  wavelengths  are  quite  close  to  each  other,  the  angu¬ 
lar  transmission  peaks  are  34°  away  from  each  other. 

Based  on  Eqs.  (5)  and  (6),  if  we  change  b  but  keep  alb 
unchanged,  the  location  of  the  peak  changes  but  the  phase 
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FIG.  7.  (Color  online)  (a)  Heights  and  (b)  full 
widths  at  half  maximum  of  the  highest  far-field 
transmission  peak  at  different  angles  generated 
by  upper  periodic  grooves  only  and  asymmetric 
upper  and  lower  periodic  grooves.  The  working 
wavelength  is  \  =  560  nm.  The  grooves  have  the 
width  Vk9=40  nm  and  depth  h=  100  nm.  The 
width  of  the  slit  is  Ws=40  nm.  ai,  a2,  £q,  and  b2 
are  optimized  for  each  beaming  angle  using  the 
algorithms  presented  in  this  paper. 
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FIG.  8.  (Color  online)  Normalized  far-field 
angular  transmission  at  two  frequencies  when  (a) 
<3j=600  nrn,  bl  =  500  nm,  a2=4 00  nm,  and  b2 
=  500  nm,  and  (b)  aj  =  540  nm,  £[=450  nm,  a2 
=  360  nm,  and  b2=A5Q  nm.  The  groove  depth  is 
h  =  500  nm. 


difference  between  the  grooves  and  the  slit  at  the  peak  angle 
don’t.  In  Fig.  8(b),  ax,  a2,  b\,  and  b2  are  10%  smaller  then  the 
corresponding  values  in  Fig.  8(a).  The  structure  still  works  as 
a  frequency  splitter  and  the  output  peaks  are  52°  away. 

In  this  paper,  we  implement  good  oblique  beaming  at  any 
angle  in  the  range  of  ±70°  using  a  metallic  sub  wavelength 
slit  with  one-side  periodic  grooves  and  better  oblique  beam¬ 
ing  at  an  angle  between  ±40°  using  a  slit  surrounded  by 


upper  and  lower  periodic  grooves.  We  also  design  a  fre¬ 
quency  splitter  using  this  kind  of  structures. 
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